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The phenomena of the microwave induced zero resistance states (MIZRS) and the microwave 
induced resistance oscillations (MIRO) were discovered in the ultraclean two-dimensional electron 
systems in 2001 - 2003 and have attracted great interest of researchers. In spite of numerous theo- 
retical efforts the true origin of these effects remains unknown so far. We show that the MIRO/ZRS 
phenomena are naturally explained by the influence of the ponderomotive forces which arise in 
1 the near-contact regions of the two-dimensional electron gas under the action of microwaves. The 

proposed analytical theory is in agreement with all experimental facts accumulated so far and pro- 
, vides a simple and self-evident explanation of the microwave frequency, polarization, magnetic field, 

■ mobility, power and temperature dependencies of the observed effects. 

Mh 

^ I PACS numbers: 73.43.Qt, 73.50.Mx, 52.35.Mw, 73.40.Cg 

I. INTRODUCTION 

- 

The discovery of the microwave induced zero resistance states (MIZRS) 1,2 in two-dimensional electron-gas (2DEG) 
systems has aroused great interest of the physical community 3 -. The MIZRS effect is observed in standard Hall-bar 
GaAs-AlGaAs quantum- well structures placed in the perpendicular magnetic field B. In the absence of microwaves 
one sees the conventional picture of the quantum Hall effect and the Shubnikov-de Haas oscillations in quantizing 
magnetic fields and the corresponding classical behavior (R xy oc B, R xx rj const) at low B; here R xy and R xx are 
the Hall and the diagonal resistance respectively. The irradiation of the samples by microwaves (with the frequency 
/ ~ 50 — 100 GHz) leads to giant oscillations of the diagonal magnetoresistance R xx at low (classical) magnetic 
fields B < 0.5 T. The effect is governed by the ratio lu/lj c of the microwave lu — 2irf to the cyclotron frequency 
u) c = eB/m*c (to* is the electron effective mass in GaAs). In the magnetic field intervals corresponding to the 
conditions k — 1/2 < uj/uj c < k, k = 1, 2, . . ., a dramatic growth of the magnetoresistance is observed (by a factor of 
i ^; 6 — 10), while at k < ui/lu c < k + 1/2 it is strongly suppressed down to "zero-resistance" states R xx s» 0, Figure []Jb). 
Exactly at uj = ku> c there is no change of the resistance, SR XX « 0. The largest amplitudes of the R xx oscillations 
are seen near the fundamental cyclotron harmonics k = 1, however the oscillation amplitudes fall down rather slowly 
with k and the oscillations remain quite visible up to k ~ 10. No apparent microwave induced changes are observed 
in the Hall resistance R xy . In the Corbino geometry very large conductance oscillations and zero conductance states 
have been also seen^. 

In the pioneering^ 2 - and the subsequent experimental papers^£~— the MIZRS effect has been observed in GaAs- 
AlGaAs quantum wells with an extremely high electron mobility /i>2x 10 7 cm 2 /Vs. A very similar effect but with 
smaller oscillation amplitudes - the microwave induced resistance oscillations (MIRO) - has been foun d 31 ' 32 two years 
earlier in the lower quality samples with the mobility /i ~ 3 x 10 6 cm 2 /Vs. In samples with fj, < 10 6 cm 2 /Vs the MIROs 
have not been observed (for the only exception see Reb^ 3 -). In such samples, instead, the photoresistance measurements 
demonstrate a completely different photoresponse^: a weak Lorenzian peak corresponding to the excitation of the 2D 
magnetoplasmon, Figure [Ha). In the first MIRO experiments^ 3 - 2 - (fi ~ 3 x 10 6 cm 2 /Vs) the magnetoplasmon 
resonance has been seen together with the w c -related resistance oscillations. 

The discovery of the MIZRS has caused an avalanche of theoretical publications^ 3 - 4 - - — . Several different scenarios 
for the explanation of the effect have been put forward. Among the proposed ideas are the so called "displacement" 
model^ ~ 36 i 39 ~ 41 i 44 , originally proposed in Refsi 68 ' 69 for a different physical situation, the microwave induced dynam- 
ical symmetry breaking 3 - 7 ., the "inelastic" model 6 - ' 47 ' 51 : 58 ' 59 , nonparabolicity effects^ 5 -, the photon assisted quantum 
tunneling 36 , phase transitions caused by electron pairing due to the exciton exchange 1 - and a quantum model that 
involves the prime number theorem 3 - 8 -. The "displacement" model suggests, for example, that microwaves assist the 
impurity and phonon scattering which leads to the spatially indirect transitions between the Landau levels N and 

N' = N + k, k = 1, 2, 3, The "inelastic" mechanism supposes that the Landau levels are substantially broadened 

due to disorder and that the microwave absorption modifies the distribution of electrons over the broadened Landau 
levels. 

The vast majority of the proposed theorie a 6 ' 34 " — has been looking for the MIZRS origin in the bulk of the 2DEG. 
They have considered an infinite 2DEG system and have examined the influence of microwaves on the bulk resistivity 
p xx of the 2D gas. We have already emphasized in Refsi 66 ' 70 that ignoring the fact that the real experimental samples 



I 



o 



- 1—^ 

X 



0.5 1 1.5 0.5 1 1.5 

(0 /CO (0 /co 

c c 

FIG. 1: This picture schematically shows the main features and the main differences of the microwave induced photoresistance 
response in (a) the moderate mobility sampled and (b) the very-high-mobility samples^ 1 ^-^. In the moderate mobility 
samples the resonance at uj c = ui^ ea is shifted from the cyclotron uj c = uj to the magnetoplasmon frequency ui 7 c es = (uJ 2 — uj 2 ) 1 ^ 2 
and the photoresistance is an even function of u c — uj^ es . In the very-high- mobility samples one observes multiple cyclotron 
resonances at u) c — ui^§ — cj/k and the photoresponse is an odd function of ui c — Here k is integer and oj p is the plasma 

frequency (J2j) . The dashed lines schematically show the dark magnetoresistance (without microwaves). 



have finite dimensions one cannot get an adequate description of the real experiments. Indeed, consider an infinite 
2DEG under the action of the external microwave field E^e~ %u)t . In the infinite sample the only relevant parameter 
which may be compared to u> is the cyclotron frequency uj c . Therefore the system response is determined by the ratio 
w/w c , for example the induced electric current 

E° 

Ut) = a xx {u>)El ex (1) 

has a resonance at the cyclotron frequency. Due to the same reason different bulk scenarios give the microwave 
induced resistivity changes determined by functions like — sin(27rw/u; c ) (e.g. in Refs t i ) or 

V- r ( E l 

> Jk const— 

k ^ Wc a 

(e.g. in RefsiS&Si^^ ^; here Jk are the Bessel functions). 

In the finite-size 2DEG samples (e.g. in the Hall bars of the width W in the x-direction) a new frequency parameter 
appears - the plasma frequency 



27rVe 2 ^ 1/2 



m*K,W 



(2) 



where n s is the 2D electron density and k is the dielectric constant of the surrounding medium. In such samples 
the internal microwave electric field E x really acting on the electrons differs from the external field E® due to the 
screening, 
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E x = i where CM « 1 - (3) 

CM ^ - ui 



is the effective dielectric function of the 2DEG stripe (see, e.g. Ref»2£). The response current ([T]) then assumes the 
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form 



3x(t) = °xx{u)E, 
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so that the cyclotron resonance is shifted to the magnetoplasmon frequency 



K+- c 2 ) 1/2 



(5) 



(the so called depolarization shift). Similarly, in all formulas for p xx suggested in the bulk-scenario papers the input 
microwave held £^ should be replaced by the screened held E x —> E®/C(u)). This would then shift all the unrelated 
resistivity oscillations and would completely destroy any seeming agreement with experiments (the depolarization shift 
is not negligible under the real experimental conditions, see a detailed comparison in Table I of Ref.— ; to the contrary, 
it is very large and has been many times observed both in the absorption^— experiments). 

Moreover, one should not ignore that (already mentioned) fact that about ten years before the hrst MIZRS exper- 
iment by Mani et ali a very similar microwave photoresistance experiment was performed in the same experimental 
group by Vasiliadou et al£. The only difference between the old (1993) and new (2002) experiments was the mobility 
of the samples (fi < 10 6 cm 2 /Vs in 1993 and /i~2x 10 cm 2 /Vs in 2002); all other parameters, such as the electron 
density, the sample dimensions, the microwave frequency, the temperature, were the same. But the results of these 
two experiments turned out to be amazingly different. Vasiliadou et al. observed a weak magnetoplasmon resonance 
([5]), Figure Ufa), which results from the screening of the external held ([3]) and is well understood. It was natural 
to expect that in the higher-mobility samples the magnetoplasmon resonance will be only slightly modihed, e.g. its 
linewidth will be reduced. In contrast, the giant w c -related oscillations have been discovered, Figure [TJb). A cor- 
rect and comprehensive theory should evidently answer the questions: Why the well-known screening effect does not 
manifest itself in the MIZRS experiments? Why the standard magnetoplasmon picture is easily observed in samples 
with relatively low mobility, while the very-high-mobility samples demonstrate a completely different response? Why 
the giant R xx oscillations and zero resistance states are not seen in low-mobility samples? These questions cannot 
be answered within the bulk models, since the depolarization shift is a consequence of the finite width of the sample. 
The main puzzle of the MIZRS phenomenon have not thus been even considered in the bulk-scenario papers. 

The bulk theories of MIZRS /MIRO have been also shown to be in disagreement with experimental facts obtained 
in the very important paper of Smet et alJ^. Being strongly related to the optical transitions between the Landau 
levels, the "bulk" scenarios are intrinsically sensitive to the sense of the circular polarization of the incident microwave 
radiation. If the "bulk" approaches were valid, the MIZRS effect, like the absorption, had to be seen at the active 
circular polarization of radiation and completely disappear at the inactive polarization. In the experiment of Smet 
et alji£ it has been found that the amplitudes of the microwave induced resistance oscillations at the right- and 
left-circularly polarized waves are essentially the same, while the absorption spectra showed a very large difference. 

The puzzle of the MIZRS/MIRO phenomena has thus been remaining unexplained so far. 

Before proceeding to the presentation of our theory let us enumerate a few further specific features of the discussed 
phenomena. 

First, the MIZRS/MIRO effects arc observed under the conditions 



which suggest that they must have a classical origin; here r is the momentum relaxation time extracted from the 
mobility, T is the temperature, Ep is the Fermi energy and N is the number of occupied Landau levels. Second, the 
MIRO are seen around u) ~ kuj c with the harmonics number k up to k ~ 10. It is known that the inter-Landau-level 
transitions N —> N' are usually strictly forbidden except for the case N N±l. This selection rule could be violated 
only if the ac electric field acting on the electrons was strongly inhomogeneous on the scale of the cyclotron radius R c . 
The value of R c under the actual experimental conditions varies from one to a few microns which is about two orders 
of magnitude smaller than both the wavelength of radiation (> 1 mm) and the sample dimensions (> 0.5 mm). This 
completely excludes the chance to bypass the selection rule N' — N = ±1. The idea that this rule could be broken 
by impurities and disorder can hardly save the situation in the ultraclean samples where disorder (if essential at all) 
should be smooth. 

In addition, it is known that in the strongly inhomogeneous external fields one should expect the Bernstein modest 
at the frequencies corresponding to the intersection of the magnetoplasmon |JSJ| with the cyclotron harmonics w = kui c . 
Usually, the Bernstein modes are very weak and seen, if at all, at k ~ 2, 3. In the MIRO-MIZRS and other similar 
absorption/magnetotransport experiments no indications on the Bernstein modes have been observed. 

Among the further puzzles of the MIRO/ZRS phenomena one should mention the observation, under certain 
conditions, of the apparent negative resistance in the finite intervals of B, Rcf. ? , and the suppression of the MIRO 



/i/r < T ~hu c <hw<s:E Fl N ~ 50 - 100, 



(6) 





FIG. 2: A schematic view of the Hall-bar sample. The circles and semicircles show electrons rotating around the cyclotron 
orbits in the bulk (6) and in the near-contact (c) regions, as well as the skipping electron orbits near the edge (e) of the sample. 

oscillations by the parallel magnetic field B\\ , ReiM-. The latter effect is observed in the quite moderate fields B\\ ~ 1 
T, when the spin effects are still not to be expected. 

Summarizing the state of the art in the MIRO/ZRS research one has to admit the presence of a very large number 
of clear experimental data, which appear to be, however, in contradiction with common sense, and the absence of 
a theory which could reasonably explain them. On the other hand, the absolute clarity and the large amplitude 
of the observed oscillations suggest that the true explanation should actually be very simple and self-evident. Such 
explanation is given in the present paper. 



II. EXPLANATION OF THE MIRO/ZRS EXPERIMENTS 

1. There exist three types of electrons in a Hall-bar sample, Figure in the bulk ('6'), near the edge (V) and 
in the near-contact regions ('c'). In the absorption experiments the main contribution to the measured signal is 
given by the bulk electrons l b\ They move in the screened electric field §3§ and absorb the microwave energy at the 
magnetoplasmon frequency ([5]) . Since the field in the bulk is weakly inhomogeneous on the cyclotron radius scale, the 
absorption spectra do not demonstrate any nonlocal resonances (Bernstein modes). 

In the transport experiments made on the samples with a moderate electron mobility /i < 10 6 cm 2 /Vs the mag- 
netoresistance response also shows the resonance at the magnetoplasmon frequency ([5]). Figure QJa). This evidently 
points to the bulk mechanism of the photoresponse in this case: electrons resonantly absorb radiation at u = w mp , 
the temperature of the electron gas increases and the bulk resistivity of the 2DEG changes. 

The dramatic change of the photoresponse in the ultraclean samples suggests that, in addition to the conventional 
bulk photoresistance mechanism, a new mechanism comes into play in such samples. This is supported by the fact that 
in the MIRO regime one can sometimes observe both the magnetoplasmon resonance and the unrelated oscillations. 
Therefore we do not attempt to find the reason of the new effect in the scattering mechanisms of the '6' electrons, 
but will search for another, bulk-unrelated contribution to the measured magnetoresistance of the sample. 

2. The origin of the second contribution to the photoresistance signal may lie near the edge or in the near-contact 
regions of the sample. The edge mechanisms of MIZRS and MIRO have been discussed in Refsi 66 i 67 . They cannot 
however explain all experimentally observed features since the edge electrons 'e', Figure [2] are practically under the 
same conditions as the '6' electrons (for example the field near the edge is weakly inhomogeneous on the cyclotron- 
radius scale, therefore it would hardly be possible to explain the higher cyclotron harmonics in the MIZRS effect). 
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FIG. 3: The electric field inside the gap between the contact wings, in the absence of the 2D electrons, at r/ c = 2ira c / 'c^fn — 100 
and W/X = 0.1 (in a typical experiment the wavelength of radiation A ~ 3 mm and W ~ 0.3 mm). 256 Fourier harmonics have 
been used, (a) and (b) - electric field is perpendicular to the boundary 2DEG - contact, (c) - electric field is parallel to the 
boundary 2DEG - contact. Inset in (a) shows the geometry of the 2D stripe between the two contact wings. 



Consider the near-contact electrons 'c'. They also feel not the external but the screened microwave electric field. In 
the very vicinity of the contacts, however, the field is screened not by the electrons of the 2DEG but by those of the 
metallic contact. Since the electron density in the metal is many orders of magnitude larger than in the 2DEG, the 
screening by the contacts in the near-contact regions is much more efficient. 

The influence of the contacts on the microwave response of a 2DEG stripe has been studied in Ref^i. In that paper 
the distribution of the electric field in the system "contact - 2DEG - contact" (see the inset to Figure [3ja)) is found 
from the solution of the Maxwell equations by expanding the fields in Fourier series over the functions ~ cos(2irxn/W). 
The contacts are described as infinitely thin 2D layers with a large (real) surface conductivity a c . Figure [3] shows the 
thus calculated electric fields E x (x,z — 0) and E y (x,z — 0) in the gap |x| < W/2 between the contacts in the absence 
of the 2D electrons (i.e. only the screening by the metallic contacts is taken into account). A very large number Nf 
of Fourier harmonics (Nf — 256 in Figure [3]) is required to get a reasonable convergency of the results for the field 
E x . Still remaining small oscillations in Figure [3ja,b) disappear when Nf is further increased. The calculations of the 
field E y parallel to the boundary 2DEG - contact converge much faster with Nf. For further details of calculations 
see RefM. 

Three important features are seen in Figure [3] 

First, the amplitude of the electric field E x (x) polarized perpendicularly to the boundary 2DEG - contact is much 
larger than that of the incident wave Eq. The field of the incident electromagnetic wave is strengthened by the metallic 
contact in the near-contact region. This is a well known effect resulting from the induced charge accumulation near 
sharp edges of metallic objects, Ref»Z&, §3. If the metallic layer i s infinitely th in and the conductivity of the metal is 
infinitely large (an ideal metal) the field E x (x) diverges as \j^J (W/2) 2 — x 2 near the edges, Ref, 7e , §3. In the real 
system with a finite a c and a finite thickness of the metallic contact in z direction the divergency is cut off but the 
near-contact field E c = E x (x = ±W/2) will still be much larger than the field of the incident wave, E c So- 
Second, the near-contact electric field E x (x) is strongly inhomogeneous on the cyclotron radius scale. In the 
experiments R c varies between 1 and 10 /xm, while the only length I relevant to the cut-off of the field divergency is 
of order of 0.1 /xm (the thickness of the AlGaAs layer on top of the 2DEG). This explains the fact that the very large 
number of the cyclotron harmonics is seen in the MIRO/ZRS experiments. Notice that the field E x (x) in the bulk, 
as well the field E y (x) in the whole sample are quasi- uniform on the scale of R c . 

Third, the near-contact electric field is strongly linearly polarized near the contact. The field E y , parallel to the 
boundary contact - 2DEG, is smaller than the incident field amplitude E . It is screened by the self- inductance effect 
caused by the tangential current j y induced in the contact by the incident radiation. In contrast, the E x component 
of the field is much larger than Eq in the near-contact region. That is, independent of the sense of the circular 
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polarization of the incident wave the near-contact electrons feel the linearly polarized ac electric field. This explains 
the inscnsitivity of the MIZRS effect to the sense of the circular polarization of microwaves observed by Smet et alJ^. 

3. Thus, near the contacts the 2D electrons move in the strong, strongly inhomogeneous and linearly polarized ac 
electric field. What physical effect can then lead to the observed photoresistance features? Remind that the MIZRS 
effect is seen in the very-high-mobility, ultraclean samples. Therefore we assume that the electron scattering is not 
important in the discussed phenomena at all and that a collisionless effect should be responsible for the MIRO/ZRS 
phenomenon. This idea is supported by comparing FiguresQIa) andQJb). The magnetoplasmon absorption resonance 
in Figure HJa) is described by an even function with respect to the resonance frequency lo = Lo mp , i.e. it is determined 
by the real part of the effective dynamic conductivity of the system (j' xx {oj). In contrast, the MIRO/ZRS resonances are 
evidently described by odd functions of lo — kuj c , Figure QJb), which strongly suggests that the considered phenomena 
are related to the imaginary part of the conductivity u' xx {lo). 

4- What could this a xx -related effect be? It is known that in the inhomogeneous oscillating electromagnetic 
field electrons experience a nonlinear, time-independent ponderomotive force, see e.g. Refs^Zr— , proportional to the 
gradient of the squared electric field. Usually the ponderomotive forces are observed in the very intense laser fields. 
To the best of our knowledge, they have never been seen in the 2DEG systems. In the absence of scattering and the 
external magnetic field the ponderomotive force acting on a charged particle is F pm (r) = — VJ7 prn (r) where 

M'l^^t)),^^,*)), (7) 

is the ponderomotive potential and (. . .)t means the averaging over time. The potential U pm {r) is proportional to 
the imaginary part of the dynamic Drude conductivity c"(w) per particle. In zero magnetic field the ponderomotive 
force always directed from the areas of the large field to the areas of the weaker field (independent of the charge of 
the particles). 

In the presence of the magnetic field B = (0,0, B), the linearly polarized electric field E(r, t) = e x E x {x) cos Lot, as 
well as a weak scattering (7 <C lo,lo c ), Eq. ([7]) is generalized as follows 

U pm (x) = ^f^(E 2 (r,;)) t 



e E x {x) I lo — lo c lo + lo c 

8m*w c I {lo — lo c ) 2 + 7 2 {lo + lo c ) 2 + 7 2 



(8) 



where 



in s e 2 10 + ij 

o~xx 0*0 = — —-, — -— -T2 2 (9) 

m* [lo + i7) z — 

is the diagonal dynamic Drude conductivity of the 2DEG [Eq. (jHJ is derived below in Section [TTTJ see also Ref. 7 -]. Two 
points are very important here. First, the absolute value of the ponderomotive force F pm cx d[E x {x)]/dx dramatically 
grows near the contacts, Figure [H[a,b). Second, in finite magnetic fields the force F pm may change its direction. If 
lo > lo c , electrons are pushed away from the high- field areas (from the contacts), while at lo < lo c they are attracted 
to such areas. The density of electrons in the near-contact region will then be reduced or increased dependent on the 
sign of lo — lo c . In other words, the microwaves form a depletion or an accumulation layer near the contacts (at lo > lo c 
and lo < lo c respectively), Figure [4] 

5. The microwave induced near-contact depletion/accumulation regions influence the experimentally measured 
magnetotransport coefficients. Consider the most interesting case when the ponderomotive forces form a depletion 
layer and microwaves suppress the resistance/conductance of the sample. In the Corbino geometry, Figure a), all 
the applied voltage then drops on the depletion regions and one observes the vanishing conductance of the Corbino 
disks. The effect is mainly determined by the ratio AT = n^/n® of the near-contact electron density n c s to the bulk 
density n s . 

In the Hall-bar geometry, Figure SJb), similarly, the formation of the depletion layers near the current contacts 
1 and 2 suppresses the measured voltage between the side contacts U xx = U34 ~ Use- To show this, consider first 
a uniform Hall-bar sample not irradiated by microwaves, Figure [SJa) . It is know n 82 ' 83 that the flowing dc current 
produces a strongly inhomogeneous distribution of the dc electrical potential <j){x, y) in the Hall bar, with power-law 
singularities at the diagonally opposite corners of the rectangular sample. For example, near the corner (a;, y) = (0, 0) 

4>{x, y) oc r v sin?7#, (10) 

where 

— «1, (11) 



contacts 



depletion/accumulation areas 



FIG. 4: The geometry of (a) the Corbino disk and (b) the Hall-bar sample under the intense microwave irradiation. The gray 
areas near the contacts show the microwave induced depletion/accumulation regions. 



(r, 9) are the polar coordinates of the point (x, y) and we assume that <r xy /a xx 3> 1. Figure EJb) qualitatively 
illustrates the behavior of the potential at the "north" and "south" sides of the Hall bar, tp(0, y) and 4>(W, y) (we have 
used cr xy /a xx = 5 in this Figure; in the experiment this ratio is much higher, e.g. cr xy /a xx ~ 250 in Refi, therefore all 
the features discussed here are much stronger in the real experiments) . The measured values of the Hall and diagonal 
voltages U xy and U xx in the absence of microwaves are also shown there. 

Now assume that the sample is irradiated and the depletion regions at y < w and y > L — w are formed near 
the contacts 1 and 2 [gray areas in Figure 03(c)]. Further, assume for simplicity that the density of electrons in the 
depletion regions and in the bulk of the sample are constant and equal n c s and rfi s respectively. The tangential electric 
current j y near the "north" and "south" sides of the rectangle must be continuous at the boundaries of the depletion 
regions, i.e. j y (0, u> — 0) = j y (0, w + 0), j y {W, L — w — 0) = j y (W, L — w + 0). In addition, the current j x should vanish 
in the corresponding boundary points, j x (0,w) — j x (W,L — w) = 0. This leads to a jump of the tangential electric 
field at the boundaries of the depletion regions, 

dyJo Pxx \dyj c n° \dy) c \dy ) c ' K ' 

As seen from Figure [5{d), the discontinuity of the field (IT2"j) leads to a reduction of the measured voltage U xx , 
proportional to the density factor J\f. The measured diagonal photoresistance R xx can then be written as a product 
of two terms, 

R xx ~ R b x Jsf. (13) 

The first factor R h xx is the photoresistance of the uniform sample (the bulk contribution). It is due to the resonant 
absorption of microwaves in the bulk and has a Lorentzian shape R xx ~ l/[(^ — u mp ) 2 + -f 2 } with the absorption 
maximum at the magnetoplasmon frequency ([5]), Figure Ufa). The second factor J\f describes the change of the 
measured photoresistance due to the near-contact microwave- induced inhomogeneity (the contact contribution). In 
the next two paragraphs we show that the density factor Af is almost always very close to unity; it may substantially 
differ from M 1 only in the very-high-mobility samples, i.e. under the conditions of the MIZRS experiments. In 
contrast, the magnetoplasmon resonance from the bulk contribution R h xx is always present in the measured signal. 
As a result, in the low-mobility samples one observes only the magnetoplasmon-resonance response-. In the MIRO 
experiments^ (relatively weak resistance oscillations, moderate-mobility samples) one sees both the magnetoplasmon 
resonance and the w c -related oscillations. Finally, in the extremely clean samples only the giant w c -oscillations are 
observed since the weak w mp resonance is hidden under the huge oscillations of R xx . This explains one of the main 
puzzles of the MIZRS/MIRO effects. 

Figure E^d) also shows that the measured variations of the Hall voltage U xy are very small as compared to U xx if 
Pxy/cxx 3> 1- This again agrees with the MIZRS experiments. 
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FIG. 5: (a) and (c): A rectangular sample with (j>(x,0) = ("west" contact) and <f)(x,L) — <f)o ("east" contact) (a) in the 
absence of microwaves and (c) irradiated by microwaves. The boundary conditions at the "north" and "south" sides of the 
rectangle are jx(0,y) = jx(W,y) =0. W is the sample width in the x-direction, L is the sample length in the y-direction; the 
arrows 5 and 6 show the position of the contacts; the arrows d show the boundaries of the depletion layers, (b) and (d): The 
distribution of the dc electric potential on the "north" and "south" sides of the rectangle (b) in the absence of microwaves and 
(d) under the microwave irradiation. U xy and U xx schematically show the measured values of the Hall and diagonal voltages. 



6. Now consider the density parameter J\f = n c s /nP s quantitatively. In the absence of microwaves the density n° in 
the uniform sample is 



m*T 



T 



(14) 



dx 



exp(x — z) 



z > 0, \z\ > 1, 
z < 0, \z\ > 1, 



(15) 



where 

is the Fermi integral in the 2D case and Co = Ep is the chemical potential (the Fermi energy). In the Fermi gas 
Co 3> T and n® = m* Xo I ' 7r ^- 2 • In the presence of the microwave induced ponderomotive potential U pm (x), Eq. ([8"]). the 
density of electrons becomes a function of x, 




n s (x) 



m*T 



T 



(16) 



The chemical potential C here may, in general, differ from Co, but, since the 2DEG is always connected to the contact 
reservoirs in the discussed experiments we will assume that C = Co = vrfi, 2 n^/m*. Now, rewrite the density (|TH)) in the 
form 



n s (x) 



Co \T 



1 



T 



= — 



Co 



Co 



T 



Co 

l-VT 2 {x)B 1 (n c ,T) 



(17) 
(18) 
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where the factor 



m*cuj 2 Ep 



is proportional to the power P = cEq /47r (per unit area) of the incident radiation, vp is the Fermi velocity, 

F{x) = ^ (20) 
is the electric field really acting on the electrons, normalized by the external field amplitude Eq, and 

RIO ^ - ^ I U) ~ ku >c " ± kuJc \ foU 

is a factor dependent on the dimensionless magnetic field O c = uj c /uj — eB/m*cw and the scattering rate T = 7/w [in 
the simplified expressions (fT5|) . ((5J) only the factor B\ is used; more general formulas with all Bk, k — 1,2,3,..., will 
appear later, see Eq. (I2"5j)]. 

7. Equation (fT5|) allows one to understand (i) why the ponderomotive effects have not been seen in the earlier 
experiments, (ii) why the MIRO/ZRS phenomenon can be observed only in the ultraclean samples, (iii) how the 
discussed effects depend on the microwave power and temperature and (iv) which factors favor and impede the 
observation of the MIZRS effect. The ponderomotive forces noticeably change the 2D electron density in the near- 
contact regions only if the correction 

%M =RF 2 (;c)Bl( o c! r) (22) 

in the square brackets in (fl~8f is not negligible as compared to unity. The basic factor in (|22|) is the power parameter 
V cx (eEoVp /uiEp) 2 , Eq. (fT9]). The factor {eE^VF /uiEp) here is the ratio of the energy which electrons get from 
the external field during one oscillation period (~ eE^vp/u) to their average energy (Ep). Usually this parameter is 
extremely small: taking for instance the typical MIRO/ZRS experimental values - P ~ 1 mW/cm 2 , / = lu/2tt ~ 100 
GHz, m* w 0.067mo and ~ 3 x 10 11 cm -2 - we get V ~ 1.3 x 10 -6 . Therefore the nonlinear electromagnetic 
phenomena can be seen only in the very intense fields and therefore the ponderomotive forces have not been observed 
in the 2DEG systems earlier. In the MIZRS experiments, however, the small parameter V in (|2"2"|) is multiplied by two 
very large factors. First, the field parameter F 2 (x) describes the giant growth of the electric field in the near-contact 
areas, Figure (3][a,b). Second, the B-dependent factor Bi(ft c ,T) becomes extremely large near the cyclotron resonance 
lj ~ uj c in the very-high-mobility samples, Figure [6] 

How the MIRO and MIZRS depend on the microwave power and temperature? As seen from Eqs. (jT5J) and (fT5j) 
there exist two different regimes. If the microwave power is not very strong, |£/ pm |/Co lj the density factor varies 
linearly with the power, 

T c = E c /Eq ^> 1. This is the regime of MIRO - the microwave induced resistance oscillations 31,32 . The temperature 
dependence of the measured signal in this regime is weak. If the microwave power is so strong that the parameter (|22[) 
exceeds unity, U pm /Co ^ 1) the density of electrons in the depletion regions becomes exponentially small (the Fermi 
gas becomes the Boltzmann one) and one gets into the regime of MIZRS - the microwave induced "zero-resistance" 
states. The P and T dependencies in the MIZRS regime arc described by the Arrhenius-type law 

/ ttp 2 PF 2 
AAccexp ^^. Bl (n c ,T) 

Such behavior of the signal as a function of power and temperature agrees with the experiments, see e.g. Refsji^^. 

Which factors favor and which impede the observation of the "zero-resistance states"? The transition from the 
MIRO to the MIZRS regime is the case when the parameter ([22]) becomes bigger than unity. Taking the maxima of 
J-(x) ~ T c = E c /E (at the contacts) and of Bi(O c ,r) (at uj — uj c = 7) we get the MIZRS observability condition in 
the form 

62 P El rf xfe)>l. (23) 
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FIG. 6: (Color online) The B-dependent factor Bi(f2 c ,r), Eq. l[2"T]). At the mobility fi = 2 x 10 7 cm 2 /Vs and the frequency 
100 GHz the value of j/uj is j/co ~ 2 x 10~ 3 (even smaller than shown in this Figure), therefore the resonance will be even 
sharper. 

The higher the radiation power P and the electron mobility fx, the easier is it to observe MIZRS. These dependencies 
have been clear from the very first MIZRS experiments. On the other hand, the high electron density n® and radiation 
frequency ui impede the observation of MIZRS. The influence of n° on the MIZRS effect has not been systematically 
studied, but the suppression of MIZRS at high frequencies has been recently reported by Studenikin et al. 21 . The 
w-dependence thus also agrees with the experimental facts. 

One more important factor influencing the MIZRS observability is the shape and the quality of the contacts. The 
ratio T c of the near-contact electric field to the external one will be larger if the contact is closer to the ideal conditions 
(the infinite conductivity <r c and the vanishing thickness in the ^-direction). There have been no systematic study of 
the role of the contacts in the microwave experiments on the 2DEG, but if the contacts, due to some reasons, turned 
out to be "more ideal" in the discussed sense, the observation of MIZRS could become possible even in samples 
with a moderate mobility. This may explain the experiment^ in which the giant magnetoresistance oscillations and 
zero-resistance states have been observed at /i < 10 6 cm 2 /Vs. 

The "negative resistance" observations of Ref.~ and the influence of the parallel magnetic fields 18 on the MIZRS 
can be also explained by the influence of contacts. In the first case one should take into account that the depletion 
regions are formed not only near the current contacts 1 and 2 but also near the side contacts 3, 4, 5 and 6, Figure 
EJb). Since the local properties of the contacts can be different, the dc field near the contacts is distorted and one 
can measure, in principle, slightly negative values of U xx . 

The second effect - the suppression of MIZRS by the parallel magnetic fields Bu ~ 1 T - seems to be completely 
unbelievable if to think about the influence of Bm on the properties of the 2DEG. This becomes, however, quite 
reasonable, if to assume that Bn modifies the contact properties. In our model the contacts are infinitely thin and 
have the B- and w-independent conductivity a c . In reality the contacts are three-dimensional and their conductivity 
tensor may be quite sensitive to the parallel magnetic fields of order of 1 T (at B ~ 1 T, m* ~ mo and / ~ 30 GHz 
the microwave and the cyclotron frequencies are equal). The suppression of MIZRS by the parallel magnetic fields 
can then be explained by the influence of Bu on the contact factor T c . 

8. So far we have discussed the microwave induced phenomena in the 2DEG using the simplified formula (|SJ for 
the ponderomotive potential. This approach is valid under the condition 



R c dE x (x) 
E x (x) dx 



< 1 



(24) 



when the inhomogeneity of the electric field is taken into account in the lowest order (the quasi- local approximation) . 
Eq. (jSJ) explains the resonance behavior of the density factor J\f near the fundamental cyclotron harmonics u> ~ uj c but 
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does not describe the resonances around u> ~ kcu c . In order to explain the microwave induced oscillations at higher 
cyclotron harmonics, one needs a more general, nonlocal theory of the pondcromotive forces. 

Such a theory is developed in Section IIIII It is shown there that the general expression for the ponderomotive 
potential has the form 



2 00 



fc=l 



where the factors Sfe(i7 c ,r) arc defined in (|2Tj) and 

1 

£fc(x) = — / E x (x + R c cos £) cos k£d£. 
t Jo 



(25) 



(26) 



Eq. (|25p contains all the cyclotron harmonics, with the amplitude of the fc-th term determined by the x-dependence 
of the electric field E x (x), Eq. ([26|l . In the special case of R c — > [the field slowly varies at the cyclotron radius scale, 
Eq. (|24|)] Eq. (|21))) gives es,(x) s» ^(x)^. Then only the first term with fc = 1 remains in the sum of Eq. (f2"5j) and 
we obtain Eq. flSJ. In the general, essentially nonlocal, case the potential (|25|) should be substituted into Eq. (JTTJ) . 
Then one gets the density factor J\f with oscillations at all cyclotron harmonics (x c < i? c ), 



T 



Co 



c 



r 



2 00 

1 ~ 1^-57- E [4-i(»c) - 4 + M] B k (n c , r) 

^° k=l 



(27) 



5. Equation (|27p gives a general formula for the near-contact contribution to the measured microwave induced 
photoresistance of the 2DEG. It depends on the behavior of the electric field near the contacts. To make our results 
more specific we consider the following model for the electric field distribution 



E x {x)=E c 



I 



I + x 



(28) 



This model describes the square-root divergency of the field near the contact which is cut off at a length I, I <gC R c 
(in a real sample I can be estimated as a distance between the contact and the 2D electron layer, i.e. I ~ 0.1 /xm). 
Calculating the coefficients efc(x) for the field (|28|) we get the density factor AT in the form (sec details in Section UTI|) 



T 



Co 



Co 



T 



k=l 



(29) 



where the functions Tk(z) are defined below in Section fill Bl [see Eq. 
five parameters: 



SJ)]. One sees that the factor Af depends on 



the dimensionless magnetic field 



eB 



mrcu) 



the dimensionless scattering rate (the inverse mobility) 



r = l = 



• the power parameter (proportional to the squared contact field electric E c ) 

2 



1 ( eE c vp 



• the dimensionless temperature T/Qq, and 

• the nonlocality parameter Vp/ul. 



12 



Figures [7] and [8] exhibit the influence of the different parameters of the problem on the density factor AT and hence, 
on the observed photoresistance R xx , Eq. (fT3"|). In Figure E{a) one sees the B dependencies of the measured signal at 
low mobilities, i.e. in the regime of MIRQ 31 i 32 . One observes how the almost constant value of the measured signal 
is transformed into the oscillating behavior. The calculated curves quite accurately reproduce those measured in the 
first MIRO experiment s 31 ! 32 . In Figure [T^b) the mobilities are higher and one gets into the "zero resistance" regime. 
A large number of higher harmonics is observed. Figure Eta) shows that the growth of the microwave power increases 
the width of the "zero resistance" regions and the amplitude of higher harmonics. The influence of the nonlocality 
parameter, Figure [SJb), is more complicated. Its reduction leads to smaller oscillation amplitudes at higher k but 
increases the width of the "zero resistance" region at k = 1. The overall agreement of the presented analytical theory 
with the experimental datai>2^i~— is evident. 

The development of the theory of MIRO/ZRS is thereby completed. 



III. NONLOCAL THEORY OF THE PONDEROMOTIVE FORCES 



A. General approach 



Consider the classical motion of a 2D electron in the uniform magnetic field B = (0,0, B) and the nonuniform 
oscillating electric field E(r, t) = e x E x (x) coswi. The equations of motion have the form 

f = v, (30) 



g 

m*v = — vxB-7m*v + F,(i,t)e„ (31) 
c 

where F x (x,t) = — eE x (x) coswi and the scattering of electrons is taken into account by the friction term —7V. The 
scattering is assumed to be small, 7/w -C 1. 

In the zeroth order in the electric field amplitude E x the equations (l3T)|) - (j3"Tj) have the solution 

^w-^r^jn, (32) 



r(°Ht)=( Xo )+R c ( Sin h!tl0> ( 33 ) 

w \yo J \- cos(uj c t + 4>) J K ' 

i.e. the particle rotates around the point (2:0,2/0) with the cyclotron frequency; here xq, yo, vq and <j> are integration 
constants and R c = vq/oj c is the cyclotron radius. In the final formulas for the 2DEG system one can estimate R c as 
Rc = vp/uj c , where vf is the Fermi velocity. Since the electric field is strongly inhomogeneous on the i? c -scale, the 
particle experiences different forces at the different parts of its trajectory and the net average force acting on it turns 
out to be nonzero. This time- independent force, which appears in the second order in the field amplitude E x , is the 
ponderomotive force F pm (xo) we are looking for. In order to find it we first search for the solution of Eqs. (|30p - (|3"Tj) 
in the first order in E x . Substituting x^(t) from (|33|) into the force F x we get 

F x = —eE x (xo + R c cos£) cos uit, (34) 

where £ = uj c t + <f> — ir/2. The function E x (xq + R c cos£) is a periodic function of £ with the period 2ir. Expanding 
it in the Fourier series 

00 

E x (x + R c cosQ = £k(xoy ki , (35) 

k— — 00 



we present F x as a sum of an infinite number of harmonics with the frequencies ±<j + kuj c . The first-order correction 
to the coordinate x(t) then reads 

00 A i(uj+kuj a )t 

x w (t)= Y p — - + ( w ->-«), (36) 
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FIG. 7: (Color online) The density factor M, Eq. (|29[) . as a function of the magnetic field at different mobilities of the 2D 
electrons: (a) the microwave induced oscillations in the moderate mobility regime (notice the vertical axis scale); (b) the 
formation of the "zero resistance" states in the very high mobility regime. Other parameters are: Vc = VT\ = 1, T/£o = 0.02, 
vf/loI — 8. 



where 



A, 



2m* 



Jk(<P-Tv/2) 



£fe(x ) 



(37) 



and the coefficients €k(x) in the Fourier expansion (|35|) are defined in Eq. (|26|) . 
In the next order in E x we obtain 



F x « -e cos uti^E x [x {0) (t)} + E' x [x m (t)]x {1) (t)) , (38) 
where E' x (x) = dE x {x) / dx. The first term in Eq. (|3"5|) is the first-order force which contains only the oscillating 
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terms. The time- aver aging of the second term should give the ponderomotive force. To calculate it we differentiate 
Eq. (|35|) with respect to x to get the Fourier expansion of the derivative E' x (x + i? c cos£), 

E' x (x + R cC o S 0= £ ^e^. (39) 
^— ' dx Q 

k— — oo 

Then we substitute (|31?|) and (j3"o]l into (j3"5)) . average the resulting expression over time and after some algebraic 
transformations finally get the ponderomotive force in the conventional form F pm (xo) = —dU pm (xo) / dxo with the 
potential ([25]) . 
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B. Ponderomotive potential for the model of Eq. (128[) 

If the near-contact microwave electric field is described by the model expression (|28j) the coefficients e k , Eq. (|26|) , 
are written as 

e k (x)=E c S k (l + j,^Y (40) 

Here the function 

1 f* cos kxdx 



1 r 
Sk(a,b) = - 

it Jo 



(41) 



S kM ^- 1 ^P^ 7 == ] ,a>\ b \. (42) 



y/a + b cos a; 

can be expressed^ in terms of the associated Legendre functions of the first kind (x) , 

(a 2 - & 2 r l/4 D , 



Substituting (HOI) into the general formula ([25]) , we get the ponderomotive potential in the form 

UpM = i^^ Rk ( 1 + V TJ Sfc( " c ' r) ' (43) 



fc=l 

J2 ^„ M c-2 



where Rk(a, b) = S^_ x (a, b) — Sj, +1 (a, 6). Evaluating the potential U pm (x) near the edge of the 2DEG we substitute 
x « i? c in Eq. (g3|) and get for t/£ m = E/ pm (x - R c ) 

^-i^E^(f)^^ r )' ( 44 ) 
fe=i v 7 

where 

T k (z) = S2_x (1 + S'+i (1 + z, z) . (45) 

The functions T k (z) are plotted in Figure [9] 

Equation (|34|l . together with ((TT]) . gives the expression ([29)) for the density factor TV. 



IV. DISCUSSION AND CONCLUSIONS 



The microwave induced giant oscillations of the magnetoresistance are thus caused by the nonlinear ponderomotive 
forces which arise in the near contact areas. The ponderomotive phenomena are well known in the plasma physics 
(they are used, for example, for the ion trapping, plasma acceleration, et cetera) but in the 2DEG systems they have 
been observed, to the best of our knowledge, for the first time in the discussed MIRO/ZRS experiments^ ' 31 ' 32 . This 
has become possible because of the very low scattering of 2D electrons in the GaAs/AlGaAs samples used in those 
works. 

The necessity to work with the ultraclean samples in the MIRO/ZRS experiments hampers the use of these effects 
in practical applications. We have seen, however, that the value of the ponderomotive forces in the considered systems 
also depends on the contact properties. In particular, the hard restriction on the mobility of the 2D electrons could 
be substantially softened in structures with the very thin (in the z-direction) contacts, in particular, in graphene 
systems^. The graphene thickness (~ lA) is about three orders of magnitude smaller than the inhomogeneity scale I 
that we have assumed above (I ~ 0.1 /xm), therefore the observation of the nonlinear plasma effects, in particular, the 
microwave induced ponderomotive forces, should be quite possible in the graphene systems. Further studies of such 
effects may lead to new interesting applications. 

So far we have discussed the MIZRS/MIRO effects only in semiconductor GaAs/AlGaAs systems. Recently a very 
similar phenomena have been also discovered in the 2DEG systems on the surface of liquid helium 26,27 , in which the 
electron mobility is also very high. These observations can be also explained by the ponderomotive forces but the 
interpretation of the effect needs some modifications since the experimental setups and parameters are quite different 
in the electrons-on-helium systems. In such systems the sample dimensions are large (of the centimeter scale) and the 
contacts are placed above the 2DEG plane (in fact, the photoresistance in the 2DEG-on-helium systems is measured 
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Z Z 

FIG. 9: (Color online) The functions T k (z) denned in Eq. (f45|) . 

by a contactless technique, see a typical experimental setup in Fig. 1 in Ref.— ). Under such conditions the field- 
amplification effect near the contacts is not the case and the contact-field enhancement factor T c is irrelevant. On 
the other hand, the 2D electrons are placed inside a cavity, in which the powerful microwave radiation produces a 
standing wave across the whole area of the system (like in the laser- field induced cold-ion traps in free space) . Under 
the action of the resulting periodic ponderomotive potential the electron density becomes inhomogeneous which leads 
to the observed microwave induced effects. 

Since in the 2DEG-on-helium systems the factor T c is not big, one could think that the observation of MIZRS 
effect would require a much higher microwave power as compared to the semiconductor systems, see the MIZRS 
observability conditions in Eq. (|23j) . This is not true, however, since the electron density in the electrons-on- helium 
systems is about four orders of magnitude lower than in semiconductors (~ 10 7 cm -2 vs 3 x 10 11 cm -2 in GaAs). 
The formation of the microwave induced electron traps should therefore be easily observable in the 2DEG systems on 
the surface of liquid helium. 
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